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Abstract 

The adsorption of functional molecules on insulator surfaces is of great importance to 
molecular electronics. Inteiplay between these molecules and insulator substrates is still not 
very clear. In this work, we present a systematical investigation of geometric and electronic 
properties of perylene-3,4,9,10-tetracarboxyUc-3,4,9,10-dianhydride (PTCDA) and Cgo on KBr(OOl) 
using density functional theory and non-contact atomic force microscopy. After energetic and 
structural details being discussed, electronic stmctures, e.g. local electronic density of states. 
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(differential) charge density, and Bader charge analysis, were inspected. It was found that elec- 
trostatics is the primary interaction mechanism for PTCDA and Ceo adsorbed on KBr, which 
can be further promoted by electronic hybridizations of non-polar molecules, e.g. Ceo, and the 
substrate. The electronic hybridization was suggested depending on the polarizability of the 
TT-system, which can be most likely further suppressed by introducing high electron affinity 
atoms, e.g. O, into the molecule. The adsorption site of a molecule on an ionic crystalline 
surface is, therefore, predominantly determined by electrostatics. Internal charge redistribu- 
tions at the molecule-insulator interfaces was found for the both systems. Since the interaction 
mechanism is predominated by electrostatics, it was thus concluded that aUcali-halides is a 
competitive candidate to be adopted to support low polarizability molecules, such as PTCDA, 
in future electronics. The influence of the found internal charge variations on electronic trans- 
port properties is still an open question, which may inspire new research fields of electronics. 

Introduction 

Using molecules as functional elements in electronic devices is a very interesting concept that 
has attracted a great deal of attention for several decades.-""- While tremendous progress has been 
realized, it is recognized that molecular electronics is a very complicated problem as not only 
isolated molecules must be understood, but also their assembly, interaction, and contact to the 
outside world, must be carefully investigated.-"- Since practically viable molecular circuitry is 
most likely to be self-assembled on solid substrates that are electrically insulated, a critical issue is 
how such a molecular layer is influenced by the underlying insulating surfaces. As demonstrated by 
many scanning tunneling microscopy (STM) studies," metal and semiconductor surfaces, most 
of which are conducting, substantially influence the electronic structure of molecules. 

Recently, ultrathin insulator films were successfully fabricated on metal surfaces,-"— which al- 
lows STM to image molecules on these ultrathin films. Non-contact atomic force microscopy (NC- 
AFM) was also demonstrated a feasible tool to investigate molecules on insulator surfaces.— 
Increasingly more experimental investigations have thus been focused on insulator surfaces or ul- 
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trathin films, e.g. KBr(001)^~— or NaCl(001).-2*i^i22i-22i However, only a few efforts have been 
made on the role of insulator surfaces to molecular overlayers, especially for electronic structures. 
Insulator surfaces were expected electronic inert to molecular overlayers, nevertheless, it has not 
yet been reported a clear physical picture of how molecule-substrate electronic hybridization is 
suppressed on insulators. In addition, the way that how molecules interact with insulator surfaces 
is also lack of full understanding. 

It is therefore the purpose of this work endeavoring to improve the understanding of these is- 
sues. In this work, we have systematically investigated the geometric and electronic structures of 
Cgo and perylene-3,4,9,10-tetracarboxylic-3,4,9,10-dianhydride (PTCDA) molecules (see Fig. [T]) 
adsorbed on KBr(OOl) using density functional theory calculations and NC-AFM. The KBr (001) 
surface is one of the most widely adopted insulator substrates for investigating molecular overlay- 
ers.—"— A polar molecule, i.e. PTCDA, and a nonpolar one, i.e. Ceo, were considered, since they 
were widely investigated and are representative molecules in molecular electronics.-"- Meanwhile, 
experimental literature of Cgo and PTCDA adsorbed on KBr(OOl) were availablei^'i^'^"-^^ for 
our results to compare with. 

In particular, DFT calculations with the dispersion correction (DFT-D)^~— were adopted to 
find the most likely adsorption configuration for each system and compared with our NC-AFM 
experiments. Adsorption energy and structural details were discussed. Given the adsorption con- 
figuration, electronic structures, e.g. local density of states (LDOS), were calculated to analysis the 
details of electronic hybridization, from which electrostatics was suggested the primary interaction 
between the molecules and the substrate. It was also found that the highest occupied molecular 
orbitals (HOMO) of Cgo electronically hybridizes with the Br p-state of the KBr substrate, form- 
ing a few new states; while there is no appreciable molecule- substrate electronic hybridization 
found in PTCDA/KBr, but a collective shift of molecular orbitals higher than HOMO, owing to a 
bending of PTCDA after adsorption. Further proofs, including differential charge density (DCD), 
real space distribution of wavefunctions, and structural distortions, were examined, all of which 
support the electrostatics. By investigating these two distinctly different molecules on a widely 
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Figure 1: Ball-stick model showing the atomic structures of PTCDA (left) and Ceo (right) 

adopted substrate KBr, our results shed considerable light on improving the knowledge of the in- 
terplay of functional molecules and insulator surfaces. Furthermore, we inferred several general 
features for the influences of insulator surfaces to molecular overlayers by comparing Ceo/KBr 
with PTCDA/KBr, which were expected to assist in choosing appropriate molecules and substrates 
for molecular nanoelectronics. 

Computational details 

Density functional theory calculations were carried out using the standard DFT-Projector-Augmented- 
Wave (DFT-PAW) method with Perdew-Burke-Emzerhof (PBE) functional— and a its revised 
form of RPBE functional for exchange-correlation energy and planewave basis sets with energy 
cutoff up to 400 eV, as implemented in the Vienna ab initio Simulation Package (VASP).— >^ Semi- 
core p electrons of K were treated as valance electrons. Five alkali-halide layers, separated by a 
vacuum slab equivalent to seven alkali-halide layer thickness (23A), were employed to model the 
surface. The two bottom layers were kept fixed during structural relaxations and all other atoms 
were fully relaxed until the net force on every relaxed atom is less than 0.01 eV/A. A p(2 x 3) and 
a p(3 X 3) supercell were employed to model the PTCDA/KBr monolayer and a p(4 x 4) one to 
approximately simulate the (8x3) monolayer— of Ceo on KBr. The surface Brillouin Zone of both 
categories of supercells were sampled by a 2 x 2 x 1 /:-point grid. 

A dispersion correction in the form of DFT-D2— was applied to PBE (PBE-D) and RPBE 
(RPBE-D) functionals, respectively. It was demonstrated that the RPBE-D method is the most 
suitable one for modeling molecule-metal interfaces among PBE, PBE-D, and RPBE-D. It is 
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expected that RPBE-D is also the most suitable one for molecule-insulator interfaces. In the rest 
of the paper, we thus focus on the results of RPBE-D while that of PEE were also reported for 
comparison. 

Results and Discussion 

Atomic Structures 
PTCDA/KBr(001) 

Recent literature indicates that the carboxylic-O-cation interaction primarily contributes to the 
molecule-substrate interaction for PTCDA adsorbed on alkali-halide surfaces.— We thus 
considered adsorption sites where the carboxylic-oxygens (denote 01) are near K cations. In terms 
of a single PTCDA adosrbed on KBr, three adsorption sites were selected among four adsorption 
sites and two adsorption orientations. They are Br-Top, Hollow, and Br-Top-R45, as shown in Fig. 
[21(a). Total energy calculations show that the Br-Top site is the most stable, with an adsorption 
energy of 1 .30 eV (0.60 eV for the PBE value), while that of the other two configurations are 0.80 
(0.20) eV for Hollow and 0.86 (0.26) eV for Br-Top-R45. 

The situation changes when PTCDA molecules form a monolayer. The smallest monolayer 
supercell that adsorption site Br-Top can form is a (3x3) one, but a smaller (2x3) one for the 
Hollow. A Smaller supercell means higher coverage and stronger intermolecular interaction, which 
may lead to configuration Br-Top, the most preferred one for single-molecule, unfavored in a 
monolayer. Two monolayer configurations were considered, i.e., the Br-Top site (the most favored 
among all the single-molecule configurations) and the Hollow site(the most compact among all 
the monolayer configurations), denoted as ML-BTop and ML-Hol , respectively (Fig. |2lb) and 
(c)). Since ML-BTop is a (3x3) supercell, its results, e.g. adsorption energy, were renormalized 
to a (2x3) supercell in order to compare with that of ML-Hol. Results of adsorption energy 
prefer the ML-Hol configuration, ascribed to a stronger intermolecular interaction. The adsorption 
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energy of ML-Hol, -1.27 eV (per molecule, the same hereinafter), is 0.34 eV more stable than 
the ML-BTop value of -0.93 eV. Those values calculated using PBE are -0.58 eV, 0.14 eV, and 
-0.44 eV, respectively. Configuration ML-Hol is therefore the theoretical suggested monolayer 
configuration. 

Figure |2] (d) shows a NC-AFM image of PTCDA/KBr, which was acquired at the edge of a 
PTCDA island. The atomically resolved KBr surface and molecularly resolved PTCDA molecules 
are also visible in that image. It clearly suggests a (2x3) superstructure where PTCDA molecules 
are along the [1 10] directions and the center of a PTCDA was assessed at a Hollow site. The shape 
of the supercell, the molecular orientation, and the adsorption site observed by the experiment are 
highly consistent with our theoretical results, which compellingly indicates configuration ML-Hol 
the most favorable structure of PTCDA/KBr(001). 

In configuration ML-Hol, the four 01 atoms of a PTCDA are on top of four K cations and the 
central phenyl ring of the molecule is at the hollow site of two Br anions and two K cations. From 
the top view of the ML-Hol, all C atoms of a PTCDA can be assigned into two categories, i.e. two 
arm-chair edges that contain 16 C atoms being on top of K cations (Over-K C atoms) and two C-C 
bonds, comprised of four C atoms, in the middle of a PTCDA where C atoms are over Br anions 
(Over-Br C atoms). Figure |2e) shows the side view of a fully relaxed PTCDA in configuration 
ML-Hol. The relaxed structure using either RPBE-D or PBE indicates a bent PTCDA, in which the 
vertical positions of the four 01 atoms are slightly lower than that of the almost planar perylene 
core, similar to PTCDA/ Ag(l 1 1).—"— As PTCDA is a symmetric molecule, the four 01 atoms 
shall be identical. When these molecules aggregate into a monolayer, a symmetry breaking was, 
however, found that the four 01 atoms diagonally split into two categories according to Tabled] 
One category is 0.1 A vertically higher than the other. We denote the higher one as OIH and 
OIL for the lower one, as show in Fig. l2lc). Three angles are available to reflect the bending 
of PTCDA, namely angles al, al, and a3. They are 174.8°, 180.2°, and 178.3°, respectively, 
suggesting a clear bending of the molecule. The larger value of angle a3, nearly 180°, stems from 
the higher vertical position of OIH than that of OIL, which is most likely ascribed to the two 
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Figure 2: Top views of the RPBE-D method optimized structures in single-molecule configurations 
(a) and in monolayer configurations (b) and (c), in which the two supercells are marked by the 
black dashed lines. The K and Br atoms are represented by larger purple and smaller brown 
spheres, respectively. The "triangle", "square", "pentagon", and "star" represent the adsorption 
sites right on top of a Br (Br-Top) and a K (K-Top) atom, in between of a K and its adjacent Br and 
at the center of the square made by two Br and two K atoms. The center of a PTCDA is indicated 
by a "-I-" in (c) and OIL, OIH, and 02 refer to lower- and higher-carboxylic O and anhydride O, 
respectively. Two blue thin dashed lines in (c) shows likely hydrogen bonds between two PTCDAs. 
A corresponding NC-AFM image to configuration ML-Hol was shown in (d). Side views of PBE 
and RPBE-D fully relaxed monolayer structures (ML-Hol) are shown in (e). 
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hydrogen bonds (bond length of 2.82 A and 2.92 A) formed between an OIH and two H atoms of 
a adjacent PTCDA, marked by the thinner blue dashed lines in Fig. I2c). 

Table [U shows the details of calculated molecule-substrate distances by RPBE-D and PBE. In 
terms of RPBE-D results, doit-K, representing the averaged distance between a OIL and the K 
cation underneath, is 3.36 A, which is 0.20A lower than that of OIH (3.56 A). Both distances 
are substantially shorter than the sum of vdW radii of O (1.52A) and K (2.75A), suggesting the 
molecule-substrate interaction shall be stronger than vdW interactions, which was further sup- 
ported by the comparison between RPBE-D and PBE. The PBE distances, in which vdW interac- 
tions were not considered, only 0.2-0.3 A longer than that of RPBE-D, still significantly shorter 
than the vdW radii sum. The qualitatively agreement between the two methods indicates that other 
interaction mechanisms, rather than the vdW interaction, are predominate for PTCDA adsorbed on 
KBr, as elucidated later. 

Table 1 : Adsorption energy (Eads) and structural details of configuration ML-Hol calculated us- 
ing PBE and RPBE-D. Angles al, a2, and a3 represent angles OIL-Ctr-OIL, 02-Ctr-02, and 
OIH-Ctr-OIH, respectively, in which the "Ctr" refers to the center of a PTCDA. Distances doii-K, 
doiH-K, and dctr-surf indicate the distances between OIL and K underneath, OIH and K under- 
neath, and from the center of a PTCDA to the substrate surface, respectively. 



ML-Hol 


Eads 


a\ 


a2 


as 


doiL-K 


doiH-K 


dctr—surf 


PBE 


-1.152 eV 


174.3° 


179.8° 


179.4° 


3.54A 


3.87A 


3.85A 


RPBE-D 


-2.532 eV 


174.8° 


180.2° 


178.3° 


3.36A 


3.56A 


3.65A 




Figure 3: Top views of the PBE and (R)PBE-D suggested most stable configurations, i.e. Cgg-DB 
(left) and CgQ-Hex (right). 
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Table 2: The PBE, RPBE-D and PBE-D results of the adsorption energies of Ceo-DB and Ceo-Hex 
and their differences. 



Adsorption Energy 




CeO-Hex 


CeO-Hex-CeO-DB 


PBE (eV) 


-0.04 


-0.02 


0.02 


RPBE-D (eV) 


-0.39 


-0.40 


-0.01 


PBE-D (eV) 


-0.56 


-0.63 


-0.07 



C6o/KBr(001) 

Totally 18 configurations of Cgo/KBr were considered based on four adsorption sites of the KBr 
surface, three contact positions of €50 , i-^- pentagon, hexagon, double-bond, and two molecular 
orientations. Two configurations, revealed by both PBE and RPBE-D, show significant stability 
than the rest, by at least 0.1 eV. In different from that of PTCDA/KBr, PBE and RPBE-D suggest 
two distinct most stable configurations of Cgo/KBr, respectively. Figure [3] left and right show the 
most stable configurations revealed by PBE and RPBE-D, respectively. In the most PBE favored 
configuration, denoted CgQ-DB, a double bond, equivalently a ~ 10° rotated pentagon, faces to the 
K cation underneath; while in Ceo-Hex, the most RPBE-D preferred configuration, a hexagon of 
the Cso contacts to the K cation underneath. 

Another method, i.e. PBE-D that the dispersion correction was applied to PBE functional, 
was adopted to clarify the issue of which method is more reliable. Table [2] shows the results of 
adsorption energy for CgQ-DB and Cso-Hex. Although PBE shows 0.02 eV more stable of C^q-DB 
than C6o-Hex, either of the other two methods indicates a more stable CgQ-Hex, 0.07 eV for PBE-D 
and 0.01 eV for RPBE-D. Since the dispersion corrected adsorption energy seems more reliable, 
it was believed that configuration Cgo-Hex is more stable than Cso-DB, which was also confirmed 
by a NC-AFM experiment. 

The contrast of NC-AFM images principally reflects of short and long range Coulomb interac- 
tions between sample and tip,— therefore, the corrugation of real-space total charge density of the 
sample can represent the height profile of NC-AFM measurements as a first order approximation. 
Figure |4] shows the the calculated total charge densities of these two configurations. Since the 
molecule may freely rotate on the substrate, the apparent height of Csq-DB shall be determined by 
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the lowest C-C bond of the top pentagon, which is 12 pm lower than the vertical height of CgQ-Hex. 
Comparing with NC-AFM data, shown in Fig. |4] middle, the height difference of 23 pm,— as well 
as the predominance of the "bright" molecules corresponding to the more stable Ceo-Hex config- 
uration, with comparatively few "dim" molecules corresponding to the lower and less energetic 
CgQ-DB, are consistent with the theoretical findings here. In configuration Ceo-Hex, the averaged 
C-K bond length of the lowest hexagon of Ceo is 3.53 A, and that of the double bond in configura- 
tion CgQ-DB is 3.28 A. Both bond lengths, clearly shorter than the sum of the vdW radii of K and 
C, also implies a stronger interaction mechanism, rather than the vdW interaction, as elucidated 
below. 

Electronic structure 

Local partial density of states (LPDOS) of PTCDA/KBr (in configuration ML-Hol) and Cgo/KBr 
(in configuration Cgo-Hex) were calculated using the RPBE-D method which gives essentially 
the same results to the PBE functional, as shown in Figs. [5]and|6l respectively. All energies were 
referenced to the vacuum level in both figures. Real space distribution of the HOMO and the lowest 
unoccupied molecular orbital (LUMO) of PTCDA and Ceo after adsorption were also discussed in 
this subsection. 




Figure 4: Isosurfaces (0.01 elh?) of the calculated total charge densities of Ceo-Hex (left) and 
Ceo-DB (right), in comparison with a NC-AFM image. The apparent height of the isosurface of 
CgQ-Hex is at least 12 pm higher than that of CgQ-DB, which is comparable to the experimental 
observed "bright" and "dim" spots (see middle). The colors of the isosurface are mapped to the 
height in the z direction (from surface towards vacuum). 
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Non-hybridized PTCDA on KBr 

Figure |5] shows the LPDOS of 01 (carboxylic O), C, and Br atoms of PTCDA/KBr before (dotted 
lines) and after (solid lines) the adsorption. It was found that the gap between HOMO-1 and 
HOMO of an adsorbed PTCDA is 0.3 eV larger than that of a bare PTCDA, but other gaps between 
adjacent upper MOs, e.g. HOMO-LUMO gap, do not change. Such a collective shift of MOs was 
ascribed to the adsorption induced structural distortion, i.e. a bent (~ 6.4° for RPBE-D) molecular 
backbone, as supported by a calculation of a bent bare PTCDA, which implies that PTCDA does 
not appear significantly influenced by the KBr substrate. In particular, the molecular states of 
PTCDA do not hybridize to substrate states in any appreciable way. The HOMO of the adsorbed 
PTCDA situates at roughly 0.3 eV lower than the upper-edge of the Br p band (~ -6.7 eV), not 
aligning to the band edge. The highest energy level of this band edge does not appreciably move, 
within several meV, before and after the adsorption. All these results suggest that the charge 
transfer between PTCDA and KBr is very small, in other words, no covalent bond was formed 
between the molecule and the substrate, which was confirmed by the differential charge density 
and Bader charge analysis discussed later. 

Hybridization between €50 and KBr 

Partial local DOSs of a C atom in the bottom layer of Cgo and one of the four Br anions under 
Cgo before and after the adsorption of Cgo on KBr were plotted in Fig. [6l It unambiguously shows 
a C-Br hybridized state, in contrast to PTCDA/KBr where no molecule-substrate electronic hy- 
bridization was found. The hybridized HOMO state, sitting at -6.45 eV, is comprised of molecular 
HOMO and a surface Br state, whereas bulk states of Br (light brown solid line) do not contribute 
to it. The adsorption of Cgo lowers the work function of bare KBr(OOl) by 0.1 eV, resulting in an 
upward shift of the upper edge of the occupied states. The HOMO-LUMO gaps between adsorbed 
(gray solid line) and bare (black dotted line) Cgo are slightly different, which is due to the upward 
shift of Cso's HOMO. These results, combined with the slightly downward shifted edge of the sur- 
face Br p band, suggest a charge transfer from KBr to Cgo, consistent with the electron-acceptor 
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Figure 5: Local partial density of states for a carboxylic oxygen (01) and an "averaged" C atom 
in PTCDA, and a Br anion under Over-Br C atoms in PTCDA/KBr. All energies are reference 
to the Vacuum Level (hereinafter). All states of the bare molecule and the bare substrate were 
plotted using dotted lines with shadows, while the states after adsorption are all in solid lines. 
Some molecular orbitals are indicated by the black arrows. The "E/r" here refers to the energy of 
the highest occupied state. 




E-Evac (eV) 

Figure 6: Local partial density of states of a C atom in the bottom hexagon layer of Cgo and a Br 
anion among the four Br anions under the molecule in Cgo/KBr. The styles of figure drawing are 
the same to those of Fig. |51 
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HOMO-C60 LUMO-C60 * ' LUMO-PTCDA 



Figure 7: Real space distribution of the wavefunctions of HOMO (a) and LUMO (b) of 
C6o/KBr(001); and HOMO (c) and LUMO (d) of PTCDA/KBr(001). 

nature of Ceo- 

Real space distribution of electronic states 

Electronic states near band gaps are of particular interest for molecular electronics. We thus plotted 
the real space distribution (RSD) of a few states originating from HOMO and LUMO of Ceo and 
PTCDA, to illustrate what was found above, i.e. electronic hybridization between Ceo and KBr 
and the electronically inert PTCDA on KBr. Figure ITJ^a) and (b) show the RSDs of the hybridized 
HOMO residing at -6.45 eV and a Ceo's LUMO locating at -4.86 eV, respectively. The -6.45 eV 
state is the hybridized state, consistent with the LPDOSs. Interestingly, the LUMO of Ceo does 
not hybridize with the substrate, hence it keeps its original shape and energy. On the other hand, 
as shown in Fig. \Tic) and (d), neither PTCDA's HOMO nor its LUMO electronically interact with 
KBr, again, consistent with the non-hybridized LPDOSs. 

Strength of hybridization 

An interesting question has thus arisen that why does Ceo hybridize with KBr but PTCDA does 
not? A key difference between PTCDA and Ceo is that in PTCDA, the O atoms have higher electron 
affinity than C atoms, which significantly increases the polarizability of the rest of the ;r-electron 
system (perylene core) of PTCDA by drawing electrons towards the ends of the molecule. The 
enhanced polarizability of ;r-electrons suppresses the ability of electronic hybridization, leading 
to the less reactive C atoms rather than O atoms for frontier MOs.— Nevertheless, O atoms in 
PTCDA still have possibility to hybridize with substrate Br, forming covalent-like bonds. However, 
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since both of Br and O are of high electron affinity, they are thus negatively charged leading 
to a repulsive electrostatic interaction rather than an attractive one, which makes the electronic 
hybridization between them unlikely. However, the negatively charged Br anions and/or O atoms 
and positively charged K cation give rise to a predominantly electrostatic interaction instead. These 
findings suggest that the polarizability of molecules largely determines the strength of electronic 
hybridization between molecules and alkali -halides. 

Interaction mechanism 

Interaction mechanism between the adsorbate(s) and the substrate is of great importance to an 
adsorption. In this subsection, two electrostatic interaction mechanisms, distinguished by whether 
it is electronic hybridization enhanced, were revealed for PTCDA/KBr and C6o/KBr, respectively, 
according to differential charge density, Bader change analysis, and structural distortions upon 
adsorption. All the calculations in this subsection were performed with the RPBE-D method if not 
specified. 

Differential charge density 

Differential charge density (DCD) is defined by poco = pTotal - P Molecule - Psubstrate- Figure [Eta) 
and (b) show the DCDs of PTCDA at the two slabs illustrated in the associated lower panels, 
respectively. A charge accumulation (warn colors) around the 01 and Over-K C atoms and a charge 
reduction (cold colors) of the 02 and Over-Br C atoms were explicitly shown in Fig. [8] (a) where 
the charge density slab is slightly below PTCDA molecules. An opposite tendency was found in 
the DCD of a slab close to the substrate (Fig. [8]^b)), i.e., a slight charge reduction of the surface K 
cations below 01 atoms and a significant charge accumulation around the Br anions under 02 and 
Over-Br C atoms. In a slab near Cgo (Fig. [8](c)), a large amount of charge accumulation was found 
around the lowest hexagon of Cgo which is over a K cation of KBr, while slight charge reductions 
were observable at the higher C atoms. In terms of the slab right above the substrate (Fig. Od)), 
a tiny charge accumulation is appreciable over the said K cation (at the center of the panel) and 
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Figure 8: Top views (upper) and side views (lower) of differential charge densities (DCD) of 
PTCDA/KBr(001) in slabs near the molecules (a) and close to the substrate (b); and those of 
C6o/KBr(001) in slabs just below the molecule (c) and just above the substrate (d). The colors are 
mapped by the DCD values, in unit of elk^, as illustrated by the color bars. 

relatively larger charge accumulations appear on the Br anions adjacent to the said K cation. 

All the DCDs share the same feature that charge accumulations and reductions are vertically su- 
perposed for both PTCDA and Ceo on KBr, which unambiguously suggests an electrostatic interac- 
tion mechanism for the molecule-substrate interaction. In terms of PTCDA/KBr, such mechanism 
was concluded the primarily interaction for the interface, since there was no appreciable electronic 
hybridization between PTCDA and KBr. Although electronic hybridization was found in Ceo/KBr, 
it induces a charge redistribution, largely enhancing the strength of electrostatic interactions. The 
electrostatic interaction is, therefore, still a very important part for the interface interaction. Bader 
charge analysis and the resulted surface structural distortions support these statements. 

Bader charge analysis 

Bader charge analysis was employed to quantitatively investigate the DCD-suggested charge redis- 
tribution and likely molecule-substrate charge transfer. Figure |9ta) illustrated the charge redistri- 
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Table 3: Charge variations of different categories of atoms for PTCDA on KBr calculated by the 
Bader charge analysis. "Layer 1" refers to the interfacial layer of KBr (hereinafter). 



Cctr Q-ed^e 01 02 H K Br 

PTCDA (e) +0.01 +0.26 +0.04 -0.01 -0.22 

Layer 1(e) ----- -i-0.02 -0.10 

Layer 2(e) ----- -i-O.OO -0.00 

Layer 3(e) ----- -i-O.OO -0.00 

Layer 4(e) ----- -i-0.02 -0.01 

Layer 5(e) ----- -i-0.02 -0.03 

Total +0.08 -0.08 



bution of PTCDA/KBr, in which electrons relocate from the H atoms (green balls) to the Over-K C 
atoms, and the middle C and 02 atoms are essentially neutral. Table [3] listed the charge variations 
of different categories of atoms that these "green balls" totally loses a charge of 0.22 e per molecule 
and a charge of 0.30 e is gained by those "red balls". An extra charge of 0.08 e was donated from 
the substrate, primarily surface Br anions, to PTCDA. 

These results of Bader charge analysis looks a bit inconsistent with the DCDs, however, it pro- 
vides more information if considering both results together. It can be inferred that the electrons 
surrounding the Over-Br C and 02 atoms were pushed up from the interstitial space between them 
and the surface Br anions to the above, due to the molecule-Br Coulomb repulsion, which forms 
electric dipoles normal to the surface, as shown in Fig. Sa). To response these dipoles, imag- 
ing dipoles were also induced at Br anions underneath. In other words, the charge redistribution 
among these atoms, i.e. Over-Br C, 02, and Br atoms, is perpendicular to the surface plane. Such 
molecule-Br Coulomb repulsion also gives rise to a in-plane charge transfer from Br anions to their 
adjacent surface K cations. Less positive K cations is against electrostatic interaction between K 
cations and Br anions. A straightforward solution is a charge transfer from the K cations to the 
PTCDA above, which strengthens both the in-plane K-Br and the normal-to-plane molecule-K 
interactions. This hypothesis was confirmed by the fact that a charge of 0.08 e was transferred 
from the substrate Br anions to a PTCDA. A in-plane charge redistribution within PTCDA was 
also found that electrons transfer from H atoms to Over-K C or 01 atoms. They are all consistent 
with the electrostatic interaction mechanism that electric dipoles at Over-Br C atoms and Br an- 
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Figure 9: Sketches of charge variations of PTCDA (a) and Ceo (b) on KBr before and after the 
adsorption. Green, gray, and red balls represent the electron-losing, neutral, and electron-gaining 
atoms after the adsorption. In term of (a), O atoms are shown by larger balls and C atoms are in 
small balls. Two electric dipole moments were marked on the lower panel by the two black arrows 
and the Br-K-PTCDA charge transfer was illustrated by the two thinner black arrows associated 
with an "e" for each of the both. Only the bottom three layers of Ceo were shown in balls in the 
upper panel of (b), while in the lower panel. Ceo was divided into seven lateral layers, from the 
first to the seventh. 

ions, the K-mediated charge transfer from Br anions to PTCDA, negative charges on Over-K C and 
01 atoms that over K cations, and additional positive charges on H atoms that close to Br anions 
undemeath and negatively charged Over-K-C and 01 atoms. 

The amount of transferred electrons from KBr to C50 is half to that of PTCDA, as shown in 
Table IH Charge redistribution was found throughout the molecule and the substrate, except the 
bottom layer of the substrate slab. The four Br anions under Ceo loses a charge of 0.04 e, while the 
all surface Br anions gain, overall, a charge of 0.01 e and a charge of 0.04 e was lost for all surface 
K cations. The most significant difference from the results of PTCDA/KBr is that a charge of up to 
~ 0. 1 e was gained by Br anions from K cations in each bulk layer of KBr, which gives rise to more 
stable bulk layers. In terms of Cgo, the first and third layers that contact to the substrate (red balls 
in Fig. Ob)) gain charges of 0.05 e and 0.04 e, respectively, while the second layer (green balls), 
in between those two layers gaining electrons, loses a charge of 0.04 e. The alternating appearance 
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Table 4: Charge variations of different categories of atoms for Cgo/KBr revealed by the Bader 
charge analysis. The "1" to "7" in column Cgo refer to the layers illustrated in Fig. ^h) and the 
"1" and "5" in column KBr represent the first and bottom layers of the KBr slab, respectively. 







Ceo (e) 




KBr (e) 






Total 
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2 3 
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5 6 7 1 


2 3 


4 
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C +0.05 


-0.04 +0.04 


+0.10 


-0.24 +0.10 +0.03 - 








+0.04 


K 






-0.04 


-0.09 -0.08 


-0.03 


0.00 


-0.24 


Br - 






- +0.01 


+0.07 +0.10 


+0.02 


0.00 


+0.20 



of electron gains and losses, also found for the fourth to sixth layers, lowers the total energy of 
Cgo- In addition, the transferred electrons to the "red" C atoms of the first layer attract with the K 
cation undemeath, further strengthening the molecule-substrate interaction. The attraction aside, 
it also reduces the molecule-substrate (Br-C) Coulomb repulsion that the electron losses of the 
second layer C atoms (0.04 e) and the adjacent Br anions (0.04 e), which considerably stabilizes 
the interface. These results explicitly suggest a charge transfer from the four Br anions, not from 
the surface K cations, to the C^q, which is consistent with the found C-Br electronic hybridization. 
All these facts manifest that electrostatic interaction, promoted by electronic hybridizations, is 
primarily responsible for the CgQ-KBr interaction, accordant with the DCD results. 

Structural distortion of the substrate 

Figure \T0\ illustrated the molecule-adsorption-induced structural distortions of KBr. In terms of 
PTCDA/KBr, the two negatively charged 01 atoms of a half PTCDA, shown in Fig. [TOta). effec- 
tively attract the two K cations undemeath, leading to the K cations being pulled out from their 
initial positions by 0.06 A along the directions indicated by the black arrows in Fig. [TOl a). Vertical 
distortions were also found for K cations and Br anions that the two K cations under 01 are verti- 
cally pulled up by 0.10 A(0.09 A for PBE), as shown in Fig. [Igb). The Br anion under the C-C 
bond appears tend to be away from the surface, since it is less negative, i.e. less attraction between 
it and the K cation undemeath after the adsorption. However, considering the electron densities 
of both the Br anion and the C atoms are rather high, a repulsion between them does exist, which 
results in the Br anion being overall pushed down by 0.05 A(0.02 A for PBE). Both are significant 
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Figure 10: Structural distortions of K cations and Br anions at the interface during the adsorption 
of PTCDA (a) and (b) and €50 (c) and (d) in top views (a) and (c) and side views (b) and (d). 
All larger balls, together with black arrows, were considered in discussion. The direction of these 
black arrows indicates the direction of movement for the considered atoms. 

vertical shifts compared to any anticipated rumpling of the KBr(OOl).— 

Similar to the PTCDA case. The repulsion between Br anions and C atoms pushes these four Br 
anion under the Cgo outwards, as indicated by the black arrows in Fig. [lOlc). While a substantial 
charge redistribution was found for Cgo/KBr, significant vertical distortions are expected for K 
cations and Br anions. Figure [TOld) shows that these four Br anions are vertically lowered by 
0.08 A (0.06 A for PBE), which is, again, ascribed to the repulsion between Br anions and C 
atoms. The more positive K cation beneath Cgo moves downwards by 0.06 A (0.09 A for PBE), 
owing to a weaker attraction between it and the surrounding Br anions whose lose a charge of 
0.04 e. All these distortions provide further proofs that electrostatics should be responsible for the 
primary interaction of this interface. 



Conclusion 



In summary, we have determined the most likely configuration for PTCDA/KBr(001), namely the 
0-over-K+ configuration, and a flat-lying hexagon and a tilted pentagon configuration for that 
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of C6o/KBr(001). On the basis of these configurations, electrostatics was revealed the primary 
interaction mechanism for PTCDA and €50 adsorbed on KBr, which can be further strengthened by 
electronic hybridizations of non-polar molecules, e.g. Cso- The electronic hybridization between 
halogens and C ;r-systems depends on the polarizability of the ;r-system, which can be most likely 
further suppressed by introducing high electron affinity atoms, e.g. O, into the ;r-system. Internal 
charge redistributions, especially vertical ones, were also found at the molecule-substrate interface 
for both PTCDA and Cgo- These charge redistributions can be regarded as a response to the periodic 
potential generated by Br anions and K cations at the KBr surface. Additional conclusions are in 
order, (a) appreciable electronic hybridization can be found for Cso/KBr; (b) the LUMO and upper 
MOs for both systems, as well as the HOMO for PTCDA/KBr, are almost uninfluenced by the 
insulating substrate; (c) the adsorption site of isolated molecules on ionic crystalline surfaces is 
predominantly determined by electrostatics. 

Due to the dominant electrostatic interaction mechanism, we conclude that alkali-halides is a 
competitive candidate to be adopted to support low polarizability molecules such as PTCDA in 
future electronics. In terms of high polarizability molecules, like Cgo, side groups or other intro- 
duced high electron affinity atoms were expected to help preventing or even tuning the hybridiza- 
tion. It would be thus interesting to examine the role of alkali-halides in other purely carbon based 
;r-systems, like graphene nano ribbons, and the role of substrate-induced intemal charge redistri- 
butions in electron transport properties in future studies . 
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